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^  This  report  contains  the  results  of  supplementary  work  done  related  to  the  reliabilltj 
analysis  of  Application  Specific  Very  Large  Scale  Integrated  (ASIC  VLSI)  CMOS 
circuits.  The  major  work  is  currently  being  carried  out  under  Task  N-9-5716.  The 
main  goal  of  both  tasks  is  to  determine  the  electromigration  susceptibility  of  VLSI 
circuits.  Electromigration  is  a  major  reliability  problem  caused  by  the  transport 
of  atoms  in  a  metal  line  due  to  the  electron  flow.  Under  persistent  current  stress, 
electromigration  can  cause  deformations  of  the  metal  lines  which  may  result  in 
shorts  or  open  circuits.  The  failure  rate  due  to  electromigration  depends  on  the 
current  density  in  the  metal  lines  and  is  usually  expressed  as  a  median-time-to- 
f allure  (MTF) . ~J 

This  work  focuses  on  the  electromigration  problem  in  the  power  and  ground  busses. 

To  estimate  the  bus  MTF,  an  estimate  of  the  current  waveform  in  each  branch  of  the 
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13.  ABSTRACT  (Continued). 

bus  is  required.  In  general,  the  MTF  is 
waveform,  and  not  simply  on  its  time-ave 


dependent  on  the  shape  of  the  current 
tage.  However,  a  very  large  number  of 


such  waveform  shapes  are  possible,  depending  on  what  inputs  are  applied  to  the 
circuit.  This  is  especially  true  for  DtvS1  circuits,  which  draw  current  only 
during  switching. The  correct  current  waveform  to  be  used  for  MTF  estimation 
is  one  that  combines,  (in  some  sense)  the  effects  of  all  possible  logic  waveforms. 
Considering  the  set  c!>f  logic  waveforms  allowed  at  the  circuit  inputs  as  a  proba¬ 
bility  space,  the  current  in  any  branch  of  the  bus  becomes  a  stochastic  process. 
The  expected  current  waveform  is  then  the  mean  waveform  (not  a  time-average)  of 
this  process.  Tnis  is ja  waveform  whose  value  at  a  given  time  is  the  weighted 
average  of  all  possible  current  values  at  that  time.  The  resulting  methodology 
is  called  a  probabilistic  simulation  of  the  circuit. 


In  this  task  it  was  proven  that  the  KTF  estimate  not  only  depends  on  the  expected 
current  waveform,  but  also  on  the  current  variance  waveform,  especially  in  regions 
where  the  current  density  is  high.  The  relationship  between  the  MTF  and  the  ex¬ 
pected  and  variance  current  waveforms  in  che  bus  has  been  established.  A  descrip¬ 
tion  of  how  an  estimate  of  the  variance  current  waveform  is  derived  at  che  terminals 
of  CMOS  gates  is  presented  and  an  explanation  of  the  approach  used  to  estimate  the 
variance  current  waveform  anywhere  within  a  bus  given  the  variance  current  waveforms 
at  the  contacts  to  thp  bus.  Implementation  issues  and  results  are  also  presented. 


UNCLASSIFIED 


EVALUATION 


This  technical  report  describes  a  technique  for  including  the  current  variance 
waveform  in  the  calculation  for  the  median-time-to-failure  (MTF)  for  metal  power 
and  ground  busses  due  to  electromigration  in  VLSI  integrated  circuits.  Previously 
described  work  details  the  use  of  probabilistic  simulations  in  determining  the 
expected  current  waveforms  in  the  branches  of  the  busses.  This  work  establishes 
the  relationship  between  the  MTF  and  the  expected  and  variance  waveforms  in  the 
bus.  Several  examples  are  presented  which  compares  the  current  waveforms 
obtained  by  using  SPICE  simulations  and  waveforms  obtained  using  probabilistic 
simulation  with  and  without  the  variance  contribution.  The  main  advantage  of  this 
approach  is  the  ability  to  handle  large  circuits  by  replacing  an  exponential  number 
of  deterministic  simulation  runs  with  a  single  probabilistic  simulation  run. 
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1.  Introduction 


Under  this  task,  we  performed  supplementary  work  related  to  the  reliability  analysis  *)f 
Application  Specific  Very  Large  Scale  integrated  (ASIC  VLSI)  CMOS  circuits.  The  major 
work  is  currently  being  carried  out  under  Task  N-9-5716.  The  main  goal  of  both  tasks  is  to 
determine  the  electromigration  susceptibility  of  VLSI  circuits.  Electromigration  is  a  major 
reliability  problem  caused  by  the  transport  of  atoms  in  a  metal  line  due  to  the  electron  flow. 
Under  persistent  current  stress,  electromigration  can  cause  deformations  of  the  metal  lines 
which  may  result  in  shorts  or  open  circuits.  The  failure  rate  due  to  electromigration  depends 
on  the  current  density  in  the  metal  lines  and  is  usually  expressed  as  a  median-time-to-failure 

(MTF). 

In  this  work  we  focus  our  attention  on  the  electromigration  problem  in  the  power  and 
ground  busses.  To  estimate  the  bus  MTF,  an  estimate  of  the  current  waveform  in  each  branch 
of  the  bus  is  required.  In  general,  the  MTF  is  dependent  on  the  shape  of  the  current  waveform, 
and  not  simply  on  its  time-average.  However,  a  very  large  number  of  such  waveform  shapes 
arc  possible,  depending  on  what  inputs  are  applied  to  the  circuit.  This  is  especially  true  for 
CMOS  circuits,  which  draw  current  only  during  switching.  The  correct  current  waveform  to 
be  used  for  MTF  estimation  is  one  that  combines  (in  some  sense)  the  effects  of  all  possible 
logic  input  waveforms.  If  we  consider  the  set  of  logic  waveforms  allowed  at  the  circuit  inputs 
as  a  probability  space,  the  current  in  any  branch  of  the  bus  becomes  a  stochastic  process.  The 
expected,  current  waveform  is  then  the  mean  waveform  (not  a  time-average)  of  this  process. 
This  is  a  waveform  whose  value  at  a  given  time  is  the  weighted  average  of  all  possible  current 
values  at  that  time.  The  resulting  methodology  is  what  we  call  a  piobabilistic  simulation 
of  the  circuit.  We  have  implemented  this  approach  in  CREST  [1]  which  uses  statistical 
information  about  the  inputs  to  directly  derive  the  expected  current  waveforms. 


In  this  task,  we  have  established  that  the  MTF  estimate  not  only  depends  on  the  expected 
current  waveform,  but  also  on  the  current  variance  xuavcformt  especially  in  regions  where  the 
current  density  is  high.  This  report  presents  our  findings  and  is  organized  as  follows.  In  the 
next  section,  we  establish  the  relationship  between  the  MTF  and  the  expected  and  variance 
current  waveforms  in  the  bus.  In  Section  3,  we  describe  how  an  estimate  of  the  variance 
current  waveform  is  derived  at  the  terminals  of  CMOS  gates.  In  Section  4,  we  explain  the 
approach  to  be  used  in  estimating  the  variance  current  waveform  anywhere  within  a  bus 
given  the  variance  current  waveforms  at  the  contacts  to  the  bus.  Implementation  issues  and 
results  arc  presented  in  Section  5. 

2.  Stochastic  Current  Waveforms  and  the  MTF 

Consider  a  metal  line  of  uniform  width  and  thickness  carrying  a  constant  current.  The 
relationship  between  the  MTF,  *50,  due  to  electromigration  in  the  metal  line  and  the  current 
density,  j,  has  been  extensively  studied,  and  shown  to  be  a  complex  nonlinear  function  [2], 
as  shown  in  Fig.  1.  We  will  consider  the  MTF  to  be  0  oc  1  //(j)  where  j  is  in  /1/cm2,  and  / 
is  a  dimensionless  nonlinear  function,  whose  plot  is  shown  in  Fig.  2,  which  was  derived  from 
Fig.  1. 

If  a  metal  line  carries  a  varying  current,  of  density  /(£),  then  the  MTF  is  £50  ex  1/Jcfl-, 
where  Jck  depends  both  on  /  and  on  the  waveform  shape  of  j[t).  It  has  been  suggested  [3] 
that,  if  the  waveform  is  periodic  with  period  T  and  consists  of  a  train  of  pulses  k  =  1 , . . . , m 
of  heights  jh  and  duration  /*,  then  : 

m  . 

«4ir  =  7pftik)-  (J) 

fcsl 

For  a  general  periodic  waveform,  we  take  the  summation  to  the  limit  and  write  : 

T 

Jcir  =  j.  /  f(W-  (2) 

0 

If  the  current  waveform  is  not  periodic,  then  better  estimates  of  Jc(f  are  obtained  by  using 
larger  values  of  T  so  that  more  features  of  the  waveform  are  included.  Therefore  one  can 
write : 

T 

Jcff  =  ^  J  f(j)dt.  (3) 

0 


2 


Now  suppose  that  the  current  waveform  is  stochastic ,  i.cM  it  is  a  stochastic  process  j (t), 
that  represents  a  family  of  deterministic  (real)  current  waveforms  _/*(£),  with  associated 
probabilities  Pkt  k=  1,...,  A\  over  the  (finite)  interval  (0,<ol-  Based  on  this  information,  we 
can  build  a  (non-stochastic)  current  waveform  j(i),  over  [0,2')  as  7’  —  sc,  that  is  indicative  of 
the  current  during  typical  operation  as  follows.  Consider  a  random  sequence  of  the  waveforms 
jjt(t),  each  being  shifted  in  time,  spanning  an  interval  of  length  <o»  and  occurring  with  its 
assigned  probability  P*,  as  shown  in  Fig.  3.  Let  nk(T)  be  the  (integer)  number  of  occurrences 
of  the  waveform  jk(t)  in  [0, 7'],  and  let  wp  =  \T/to\.  If  Jk,  =  l,...,Af  are  defined  as  follows  : 

<0 

4=jj//C#»R  ('0 

°  0 

then  : 


By  the  law  of  large  numbers  (4),  liin„r_.M(7ij.(21)/nT]  =  Pt,  which  leads  to  : 

Jctr = t  [£  J  nh)dtj  pk = i  / ! \t  mPt]  dt  (o) 

and,  finally  : 

*D 

J*  =  5  /  £|/(j)]<«  (?) 

0  0 

where  JE( )  denotes  the  expected  value  operator.  This  is  an  important  result;  it  says  that 
the  MTF  due  to  a  stochastic  cuirent  depends  only  on  the  expected  wavefonn  of  a  nonlinear 
function  of  the  current. 

Since  f  is  nonlinear,  J5[/(j)]  is  not  easy  to  evaluate.  At  low  current  values,  where  /  is 
linear  (Fig.  2),  i3[/(j)]  =  /(JSfj)).  If  this  is  substituted  in  (7)  and  compared  with  (2)  it  shows 
that,  when  /  is  linear,  the  expected  current  waveform  JS[j]  derived  in  (1)  may  itself  be  used  as 
the  current  waveform  j(t)  in  (2)  for  MTF  estimation.  This  establishes  the  importance  of  the 
expected  current  waveform  for  electromigration  failure  analysis.  In  general,  /  is  nonlinear, 
and  a  generalized  approach  will  be  developed  below. 


At  any  time  i,  the  process  j(f)  can  be  thought  of  as  a  random  variable  j  with  mean 
j/j  =  £(j),  and  variance  aj  =  iS[(j — In  general,  the;)'*  moment  of  j  is  ftJ})  =  £?[(j— 

To  estimate  the  mean  of  /(j),  £(/(j)),  we  use  a  Taylor  series  expansion  of  /,  which  leads  to  : 

£|/(i)| »  fi’ij) + /'%)$  +  •  •  • + (S) 

-  r* 

It  is  evident  that,  when  /  is  linear,  (8)  reduces  to  : 

£(/()))  =  /(BUI),  (9) 

as  observed  above.  Hence  using  the  expected  current  waveform  as  an  actual  current  waveform 
for  MTF  estimation  based  on  (2)  amounts  to  making  a  first-order  approximation  in  (8). 
Naturally,  higher  order  approximations  would  lead  to  better  results.  In  particular,  if  /  is 
approximated  by  a  quadratic  in  the  neighborhood  of  i/y,  then  : 

mawfa) +/%)$•  (i») 

This  second-order  approximation  becomes  exact  if  f(j)  is  represented  by  the  straight  lines 
corresponding  to  j 1  and  j2  in  Fig.  2.  It  is  more  accurate  than  (9)  since  it  covers  a  wider  range 
of  currents.  As  a  result,  equations  (10)  and  (7)  offer  a  new,  more  accurate  technique  for 
computing  the  MTF.  In  order  to  make  use  of  this  technique,  we  need  to  derive  the  variance  of 
the  current  waveform  in  addition  to  its  expected  value.  As  pointed  out  in  the  introduction, 
the  estimation  of  the  expected  current  waveform  has  already  been  described  in  our  previous 
work  [I];  the  next  section  will  discuss  the  derivation  of  the  variance. 

3.  Derivation  of  the  Gate  Variance  Waveforms 

Figure  4  shows  a  generic  CMOS  gate  structure  which  we  will  use  in  our  derivation.  The 
variance  waveforms  for  the  gate  total  and  output  currents  will  be  modeled  by  triangular 
pulses  K[t<0<(t)]  and  F[i(0),  respectively,  with  peak  values  of  V[It0t]  and  V[/J.  If  an  event 
occurs  at  the  gate  input  at  time  t,  then  we  denote  by  t~  and  t+  the  instances  of  time 
immediately  before  and  after  the  event,  respectively.  Focusing  for  now  on  the  output  current 
pulse,  its  variance  waveform  starts  with  a  peak  of  V[f)  =  V'[t(t+)]  at  lime  t  and  decays  linearly 
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to  zero  at  time  d  +  t.  Since  V'[/)  =  E(/2]  -  E(/)2  (4),  and  since  CREST  already  derives  the 
expected  pulse  peak  (El/)),  we  will  concentrate  here  on  the  derivation  of  E[/2]. 

Let  ip  =  ip\  +  ip 2  and  i„  =  t„,  +  i„ 2-  It  is  easy  to  verify  that  ipl  =  ip  x  C„/(C,,  +  C„), 
and  —  in  x  Op /{Op  4*  Cm)-  lliercfore  ! 


SpW)  =  JB[i?W) 


£$(0) 


The  term  containing  £(ip(l)>n(0)  's  omiUal  (it  is  zero)  since  at  least  one  of  the  charging 
currents  is  zero  at  any  given  time.  In  particular,  the  value  at  the  peak  is : 


E[m)  =  £(/;«)  a 


+  £(/2(t)|  x 


The  values  of  £(/£]  and  13(7; |)  arc  derived  as  follows.  For  E[/,7),  consider  the  p*part  of 
the  gate,  and  let  every  transistor  T*  be  represented  by  a  switch  of  on-conductance  <70„,jt  (5). 
Based  on  this  switch-network  model  of  the  p-block,  let  GP(t)  be  the  random  conductance 
between  the  output  node  and  V^a.  Gp  is  a  function  of  the  individual  transistor  random 
conductances  <7*,  where  <7*  is  0  if  the  transistor  is  off  and  g0Ui j.  if  it  is  on.  If  an  event  occurs 
at  the  gate  at  time  d,  then  the  value  of  C?r(d+)  and  the  previous  state  of  the  output  node, 
V<>(i" ),  will  determine  Ip.  Formally,  we  have  E[/£]  =  E((V^j  -  V0(d~))2  x  G2(d+)j,  whicli 
becomes  : 

E[/2)  -  V£  x  E[G,2(d+)  1  Gp(l~)  =  0)  x  />(G„(i-)  =  0) 

where  P{A)  is  the  probability  of  the  event  A,  and  E[A  \  B\  denotes  the  conditional  expected 
value  (4)  of  A  given  B.  The  formula  is  correct  because  if  Gp(l~)  —  0  (jfc  0)  then  K>(d“)  =  0 
(K/j)-  Similarly  for  the  n-part  of  the  gate,  we  get : 


E{/2]  =  Vl  x  E[Gl{i+)  |  G„(l-)  =  0)  x  P(G„{t -)  =  0) 


To  derive  the  conditional  expectations,  consider  a  graph  representation  of  the  p-block 
(or  n-block),  where  every  edge  in  the  graph  is  labeled  with  E[gj:(t+)  1  Gp[i~ )  —  0), 
£[^.(i+)  |  GP{t~ )  =  0},  and  the  gate  node  probabilities  of  its  corresponding  transistors. 
The  details  of  how  these  quantities  can  be  derived  for  every  transistor  can  be  found  in  (5). 
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Then  perform  a  graph  reduction  operation,  which,  simply  stated,  involves  a  number  of  se* 
ries/parallcl  combinations  and  node  eliminations  that  reduce  the  graph  to  a  single  edge, 
whose  labels  are  the  required  statistics  £(C?"(£+)  j  GP(t~)  -  0]  and  E[GP[t+)  |  C7,,(t~)  »  0]. 
Similarly  foe  the  n*block. 

Having  found  the  peak  V[/J  =  J G(/2]  -  E[I ]2  for  the  output  current,  the  time  span  r 
will  be  found  by  first  solving  for  the  area  under  the  V[»(0)  pulse.  Notice  that,  if  i{t)  is  a 
triangular  pulse  of  height  I  and  area  qt  then  : 

In  this  case,  q  is  equal  to  the  charge  delivered  to  (or  from)  the  output  node  capacitors.  From 
this  it  follows  that : 

/”  =  Ill'll,  ami  jT  £[;(t)]2<l!  = 

The  second  equation  follows  since  E[i(l) ]  is  a  triangular  pulse  of  height  E[l\  and  area  E[q], 
Therefore,  tiro  variance  pulse  has  an  area  : 

= jf  V(, •(«)]<«  =  §(£[/«)  -  £(/!%!). 

The  value  of  E[Iq)  can  be  written  as  : 


£(Jr/]  =  *1-  /„i)  x  {qp]  + 


where  /,, j  (/„i)  is  the  peak  of  ipi(t)  (i„i(t)),  and  qpl  (qn j)  is  the  charge  delivered  ly  ip\(t) 
(lni(0}*  When  the  p-block  is  conducting,  the  n-block  is  not  conduling,  and  >V,(t)^0, 

-  0,  qP[  «  VjjCn,  and  </„,  »0;  while  when  the  n-block  u  conducting,  the  p-block  is  not,  and 
tiiiWirO.  iPl( t)  —  0,  9„,  w  VjdCPt  and  qPl  w  0.  It  follows  then  that: 


T!io  time  span  of  the  gate  output  current  variance  pulse  is,  therefore 

j/,£[/e1-g|/|S[.)l\ 

n  m  y 


G 


If  i(0((l)  is  me  total  gate  current,  then  : 


jT  I'M*  =  |  (£|/, -  £(/,.,]  EM)  . 


Unfortunately,  iT (/foi7rof)  dues  not  have  a  simple  expression  as  was  found  for  £[/f/]  above.  We 
have  chosen  to  use  a  conservative  estimate  based  on  the  following  assumption  .  whenever 
a  node  in  the  p-block  (n-block)  is  charged  to  Vdd  (V'Ja),  then  every  other  node  in  the  p- 
block  (n-block}  is  also  charged  to  I'jj  (V^).  This  assumption  is  true  for  simple  gates,  and 
overestimates  the  current  charge  product  in  more  complex  cases.  Based  on  this  assumption, 
one  can  show  (5)  that : 


pit  -  I  ^lOp.tct]  QpCn^U))]  ,  J'AQi^tot]  QnCttB[Iu) 

[htqtot J~  "^+er+  Ei<]n\  Tvfcr1 

where  E[qp ]  and  £%„)  arc  available  as  equations  (9)  k  (10)  in  (6),  and  £[<7a,<ot]  ‘ire, 

respectively,  the  first  and  second  summations  in  equation  (7)  in  [6],  and 


Q,.=  E  Qn  -  £ 

iiP  block  k-N  block 


As  was  assumed  for  the  expected  current  pulse  [5,  6]  wc  let  the  time  span  of  the  gate 
total  current  variance  pulse  be  equal  to  that  derived  for  the  gate  output  current,  therefore  : 

i/t  r  y_[E\I,ot<ltot)-E[Itot)E[qiot)' 

1 1  w^Wm 


x  V|/|. 


4.  Estimating  the  Variance  Current  Waveforms  in  the  Bus 

Since  the  current  density  j(d)  in  any  branch  of  the  power  or  ground  bus  is  directly 
proportional  to  the  current  i(t)  in  that  branch,  then,  to  simplify  the  presentation,  we  will 
discuss  the  derivation  of  oj-  rather  than  a-.  The  variance  of  the  current  waveform  is  a  time 
function  <7*(0  which  we  will  refer  to  as  the  vaviancc  wavefonn.  Furthermore,  wc  will  discuss 
the  power  bus  only  since  the  ground  bus  analysis  is  similar. 

The  current  in  a  branch  of  the  bus,  i(t),  is  a  function  of  the  currents  being  drawn  off  the 
bus  contacts,  ij(0>i  =  l,...,n.  Each  of  these  is,  in-turn,  simply  the  sum  of  the  individual 
gate  currents  tied  to  each  contact : 

ij(0  =  iji(0  +  — Mjk(t).  (11) 
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Tims,  in  she  framework  of  our  probabilistic  simulation  technique,  tlu*  process  of  deriving  the 
variance  waveforms  consists  of  three  steps  : 

•1*  Using  the  statistics  of  the  signals  at  the  inputs  to  each  logic  gate,  derive  the  variance 
waveform  for  its  current. 

-2-  Combine  these  variance  waveforms  at  each  contact  point  to  derive  the  variances  of  the 
contact  currents. 

-3-  Using  the  bus  topology,  and  the  variances  of  the  contact  currents,  derive  the  variances 
of  the  bus  branch  currents. 

Step  1  has  been  described  in  Section  3  above.  The  other  two  steps  will  be  described  below. 

The  critical  issue  is  the  correlation  between  the  different  current  waveforms.  Since  such 
correlation  is  too  expensive  to  derive  for  VLSI  circuits,  wc  will  occasionally  be  making  con* 
servativc  approximations  to  simplify  the  problem.  Our  experience  with  the  probabilistic 
simulation  approach  suggests  that  neglecting  the  correlation  between  different  current  wave¬ 
forms  gives  good  results  in  most  cases. 

Based  on  this  premise,  we  assume  that  the  gate  currents  tied  to  the  same  contact  are 
unconvktcd.  This  immediately  provides  a  simple  solution  for  step  2,  using  (11),  ns  follows  : 

^(0=%(0  +  -  +  ^t(0*  C12) 

The  remainder  of  this  section  will  be  devoted  to  the  more  difficult  task  of  solving  step  3,  i.e., 
deriving  the  bus  current  variance  waveforms  from  those  of  the  contact  currents. 

The  metal  bus  can  be  modeled  as  a  multi-input  multi-output,  causal,  linear,  time- 
invariant,  (LTI)  system  with  causal  inputs  Xj  and  outputs  yi.  The  inputs  Xj (t),j  =  1 ,...,« 
represent  the  contact  currents,  and  carry  the  stochastic  processes  ij [t)  of  known  variance 
waveforms  or?  (t).  The  outputs  yi(£)>  t  =  l,...,m  represent  the  l  us  branch  currents  at  which 
the  variance  waveforms,  cj.(i),  arc  required.  Let  be  the  impulse  response  function 
relating  yt (L)  to  Xj [t)  : 

y»(0  =  ]£%(*)  *xj(0>  i  =  (13) 

j=i 

where  denotes  the  convolution  operation. 
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It  is  well  known  (see  (*1),  page  209)  that  the  variances  of  the  system  inputs  are  not  enough 
to  derive  the  variances  of  its  outputs.  The  auto-conxlntion  of  each  input,  (£j,  £3)  = 
/?(xj(ti)xj(t2))>  is  also  required.  Since  the  input  processes  are  not  undc-scme  stationary  (<l), 
an  exact  analytical  solution  can  be  quite  complex,  even  if  the  autocorrelation  were  known. 
Therefore,  as  is  often  necessary,  we  will  make  certain  simplifying  assumptions  about  the 
structure  of  ftzjxj* 

Wc  will  assume  that  the  correlation  between  Xj(£)  and  xj (t  +  r)  goes  to  zero  as  r  —  00. 
In  terms  of  the  auto-covariance,  Cx}xj{t\yh)  =  //^(f^fa)  -  Ibis  will  be 

formulated  as  : 

C*^(£i,fj)  =  0x^(£i),  and  OxjXj(thh) »  0  for  |£j  —  £2!  £  J «  (1J1) 

where  T  is  a  (typically  small)  time  interval. 

Consider  the  discrete  lime  system  obtained  by  sampling,  with  period  T,  the  continuous 
time  system  defined  by  (13).  If  Xj(&]  =  Xj(fc!F)  are  the  discrete  processes  at  the  inputs,  and 
y\[k)  =  yi(kT)  are  the  discrete  output  processes,  then  : 

yi[*]  =  E  j  *  *#)»  *  -  1 . m  ( 1 5) 

i=i 

where  h\^[k]  is  the  discrete  impulse  response  function  relating  yi[k)  to  Xj [&).  As  shown 
below,  the  discretized  output  variance  waveforms  can  be  derived  irrespective  of  the  shape 
of  Cxjx^tijto)  for  |£j  -  £3!  <  T.  The  continuous  variance  waveforms  can  then  be  obtained 
by  interpolation.  Strictly  speaking,  therefore,  the  sampling  period  T  should  be  small  ;  1 JT 
should  be  larger  than  the  largest  frequency  component  of  the  inputs.  However,  since  fine 
waveform  details  are  not  of  paramount  importance  in  this  work,  wc  need  only  restrict  T  to 
be  small  enough  so  that  waveform  features  in  so  small  an  interval  arc  inconsequential. 

To  simplify  the  notation,  define  yy[fc]  =  h[j\k]  *  Xj[fc].  Furthermore,  as  pointed  out 
above,  we  will  neglect  the  correlation  between  the  contact  currents.  Hence  the  Xj  inputs  are 
uncorrelated,  and  : 

°%  [*]  =  E  aij  (*1»  t  =  1 , . . . ,  m. 
i=i 
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(16) 


We  have  tints  reduced  the  problem  to  analyzing  a  single-input  single-output  discrete  ITl 
system  : 

yu[*M?’W  •>#)=  E  *!>!># -4  on 


<=Q 


Let  xj [k]  =  xj[k]  -  i )Xj[k]  and  y#)  =  y#)  -  »/*,(*)•  Then  [fc]  =  £[y#]2)  mul 
y#)  =  h^lk)  *  Xj(t'),  hence  : 

<(<••] =  t  *!?h)  t  (W) 

L  '<=o  *  J  <1=0  <5=0 

Furthermore,  it  is  easy  to  see  that  i2(xj[t'i]xj(t-2)J  =  C^xj^hfcs),  which,  using  (14),  gives  : 

■ .ai,w.  os) 

<=0 

And,  finally,  the  variance  waveforms  for  the  system  outputs  arc,  using  (16)  : 

i  =  l . m.  (20) 

j- * 


In  other  words,  the  variances  of  the  system  outputs  (bus  branch  currents)  can  be  obtained 
from  the  convolution  of  the  variances  of  its  inputs  (contact  currents)  with  the  squares  of  its 
discrete  impulse  response  functions.  This  discrete  convolution  can  be  easily  performed  once 
the  discrete  impulse  response  functions  are  found.  Of  course  the  summation  need  not  be 
taken  to  infinity,  and  may  be  conveniently  truncated  after  is  less  than  some  small 

value.  To  obtain  the  discrete  impulse  response  functions,  note  that  if  a  unit-step  input 
current  is  applied  at  contact  j,  with  all  other  contact  currents  held  at  zero,  and  if  the 
resulting  outputs  y,(t)  arc  monitored,  then  : 

Jfp 

M?  W  =  y,-(*T)  -  »,-((*  - 1)3")  =  Mr)*,  ;  =  1 . (21 ) 

This  suggests  two  methods  for  deriving  The  first  uses  a  simulation  program  such  as 

SPICE  to  simulate  the  bus  with  unit-step  input  currents  applied  at  each  contact  (one  at  a 
time),  while  monitoring  the  bus  branch  currents.  This  gives  the  nm  functions  h[j  *(&]  using 
(21).  Another  (approximate)  method  would  be  to  make  use  of  the  second  equality  in  (21)  . 
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if  the  continuous  impulse  response  functions  Are  Approximated  using  some  RC  time-constant 
Analysis  of  the  bus,  then  the  discrete  impulse  response  functions  can  be  obtained  from  them. 

For  very  large  chips,  it  may  be  prohibitively  expensive  to  perform  the  required  convolu¬ 
tions.  One  can  simplify  the  calculations  by  making  an  additional  assumption  as  follows.  If 
the  bus  is  known  to  be  “fast",  i.e.,  if  dies  down  faster  than  changes  in  cr£  (Jr),  then 
(19)  reduces  to  : 

(22) 

*&0 

So  the  convolutions  in  (20)  can  be  replaced  by  simple  multiplications,  and  the  constants 
can  kc  derived  in  a  pre-processing  step  from  the  impulse  response  functions 
and  stored  in  a  single  m  x  n  constant  matrix. 

If  the  chip  is  too  big  to  even  derive  then  one  further  simplification  can  be  made 

as  follows.  If  /t^(Jt)  dies  down  faster  than  changes  in  Xj[Jfc]  then  (17)  reduces  to  y*j(A:)  = 

xjIMds.oMJ’m-  ,uu|  s°  ; 

(23) 

The  constants  (E«s0^if  M)2  can  vcry  easily  obtained  as  follows.  Note  that  *1?M 
is  the  steady  state  current  in  branch  i  in  response  to  a  unit-step  input  current  at  contact  j, 
with  all  other  contact  currents  held  at  zero.  If  the  bus  is  modeled  as  a  resistive  network, 
then  the  steady  state  node  voltages  in  response  to  such  inputs  arc  the  entries  of  the  driving 
point  impedance  matrix.  So  if  the  node-admittance  matrix  is  built  by  simple  inspection  of 
the  bus  and  then  inverted  to  produce  the  driving  point  impedance  matrix,  the  steady  state 
currents  are  immediately  available. 

5.  Implementation  and  Results 

The  variance  calculation  technique  outlined  above  has  been  implemented  in  CREST.  YVe 
present  below  the  results  of  CREST  runs  on  a  variety  of  circuits,  showing  both  waveform 
comparisons  and  timing  performance. 

We  start  out  with  a  simple  example,  a  2-input  CMOS  NAND  gate.  The  variance  wave¬ 
form  comparison  between  CREST  and  SPICE  is  shown  in  Fig.  5.  The  SPICE  waveform  is 
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derived  by  running  SPICE  on  the  NANI)  gate  for  all  possible  logical  transitions  ai  its  in 
puls,  deriving  the  expected  current  waveform  by  doing  a  lime-point  o-eraging  of  the  results, 
and  then  using  that  to  find  the  variance  as  the  lime-point  average  of  (t  -  J£(»))2.  Since  the 
object  of  this  research  is  to  handle  very  large  chips,  ami  since  elect,  .migration  models  fur 
ac  waveforms  arc  still  controversial,  it  makes  little  sense  to  shoot  for  pc  feci  eeuracy  in  the 
current  waveforms.  It  is  more  important  to  be  able  to  derive  in  a  very  abort  lime  a  .  acefurm 
that  matches  the  peak  and  general  shape  of  the  SPICE  waveform. 

Another  single-gate  comparison  is  shown  in  Fig.  6  for  a  CMOS  complex  gale.  The 
comparisons  for  two  larger  circuits  are  shown  in  Fig.  7  (for  an  XOR  circuit)  and  Fig.  S  (fur 
a  5-l-MOSFET  2-bit  ripple  adder  circuit). 

Our  final  example  is  a  much  bigger  and  more  complex  circuit;  it  is  a  64S-MOSFET  Tbit 
parallel  multiplier.  This  circuit  is  too  big  to  make  the  2s  required  SPICE  simulations.  \\V 
will,  therefore,  show  two  different  CREST  runs  to  demonstrate  that  the  variance  computation 
works  well  even  when  the  heuristics  introduced  for  handling  large  circuits  in  (1)  are  used. 
In  Fig.  9  we  compare  a  full  accuracy  CREST  run  and  a  heuristic  CREST  run  in  which 
all  internal  nodes  of  the  multiplier  were  assumed  independent.  The  excellent  agreement 
reaffirms  the  conclusion  made  in  (l]  that  as  a  circuit  becomes  larger,  the  correlation  between 
its  internal  nodes  may  be  safely  neglected. 

YVc  next  examine  the  importance  of  the  variance  waveform  for  MTF  estimation.  The 

expected  and  variance  waveforms  combine  to  provide  a  dc  current  density  value  tu 

represent  the  ac  current  waveform  for  MTF  estimation  as  given  in  (7)  and  (10).  Equation 

(7)  is  first  evaluated  using  only  the  first  term  of  (10),  i.c.,  using  only  the  expected  waveform, 

to  give  Jc f[|.  Equation  (7)  is  evaluated  again  using  both  terms  in  (10),  i.c.,  using  both 

expected  and  variance  waveforms,  to  give  The  percentage  contribution  of  the  variance 

waveform  to  J.rr  is  then  measured  as  x  100.  The  results  are  tabulated  in  the  fifth 

c,s  Jm 

column  of  Table  l  for  a  number  of  test  examples.  The  results  clearly  establish  the  importance 
of  the  variance  waveform  in  addition  to  the  expected  waveform. 

Finally,  we  illustrate  the  speed  performance  of  CREST  with  the  variance  estimation 
built  in.  Table  1  shows  also  the  speed  comparisons  between  CREST  and  SPICE  for  all  the 
examples  presented  above.  The  speedup  becomes  much  better  for  larger  circuits  (1529X  fui 
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the  adder  and  1 1595X  for  the  multiplier).  In  fact,  the  speedup  should  grow  exponentially, 
because  an  exponential  number  of  deterministic  simulation  runs  are  replaced  by  a  single 
probabilistic  simulation  run.  We  point  out  the  case  of  the  multiplier  circuit  (the  largest 
circuit  in  the  table)  with  the  heuristic  CREST  run  (last  row  in  Table  1);  considering  the 
excellent  waveform  comparison  in  Fig.  9  along  with  the  dramatic  speedup  of  11595X  in 
Table  1,  this  establishes  the  feasibility  of  solving  large  VLSI  chips. 


Table  1.  Execution  time  comparisons.  Time  is 
in  CPU  seconds  on  a  VAX-11/7S0;  size  refers  to 
the  number  of  transistors. 


Circuit 

Size 

SPICE 

CREST 

Variance 

Contribution 

Speedup 

Nand 

4 

41.75 

0.90 

110% 

46X 

Complex 

G 

244.15 

1.09 

185% 

224  X 

Xor 

IG 

45G.00 

3.13 

236% 

14  GX 

Adder 

54 

32620.42 

21.33 

107% 

1529X 

Multiplier 

CIS 

G97530.SS* 

1S71.99 

219% 

373X 

Multiplier 

CIS 

697530.SS* 

60.161 

200% 

11595X 

‘Estimated  (2*  times  the  cost  of  a  typical  logical  SPICE  run). 
I  Heuristic  CUES'!'  run,  all  others  are  full  accuracy  CREST. 


6.  Summary  and  Conclusions 

We  have  discussed  the  problem  of  estimating  the  median  time-to-failure  (MTF)  due  to 
clcctromigration  in  the  power  and  ground  busses  of  CMOS  VLSI  circuits.  In  this  task  we 
have  verified  that  including  the  variance  wavefonn  of  the  current,  in  addition  to  the  expected 
waveform  derived  in  (1),  further  improves  the  accuracy  of  MTF  estimation.  This  was  done  by 
showing  that  the  variance  contribution  to  the  MTF  estimate  can  be  in  the  range  of  100%  to 
200%  relative  to  that  of  the  expected  current  waveform.  We  have  described  a  novel  technique 
for  deriving  the  variance  waveform,  and  its  implementation  in  the  probabilistic  simulatoi 
CREST.  The  results  of  several  CREST  runs  have  been  presented,  and  they  show  good 
waveform  agreement  with  SPICE,  as  well  as  excellent  speedups  over  traditional  approaches 
-  a  speedup  of  over  11000X  was  demonstrated  on  a  G4S- transistor  circuit. 

This  work  proves  that  the  expected  and  variance  waveforms  of  the  stochastic  current 
model  are  :  (1)  essential  to  derive  an  accurate  MTF  value,  and  (2)  can  be  efficiently  derived 
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using  the  probabilistic  simulation  approach.  The  main  advantage  of  this  approach  is  the 
ability  to  handle  large  circuits  by  replacing  an  exponential  number  of  deterministic  simulatum 
runs  with  a  single  probabilistic  simulation  run.  Without  such  a  technique,  analyzing  the 
reliability  of  large  circuits  would  seem  to  be  an  impossibility. 
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Figure  l:  The  dependence  of  MTF  on  current 
density,  reproduced  for  convenience  from  (2).  The 
dashed  lines  show  the  results  of  the  approximation 
<50  oci"1*  for  n  bs  1,  3/2,  and  2. 


Figure  2:  A  plot  of  f[j),  obtained  from  Fig.  1  by 
inverting  and  appropriately  scaling  the  ordinate 
axis. 
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Current 


Figure  3:  A  (non-stochastic)  current  waveform, 
j(t),  Uuilt  ns  a  sequence  of  the  waveforms 
each  umtring  with  its  assignee]  jirubnbmty  / 


Figure  4:  A  generic  t'MOS  gate  structure 


Ar**HCc  rciu*t4  ^0r  x  3’topui  complex 
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figure  3:  Variance  results  fur  a  Sl-MOSFET  2- 
lut  ripple  adder  CMOS  circuit. 


Figure  9:  Variance  results  for  a  UlS-MOSFET 
*Uuit  parallel  multiplier  CMOS  circuit. 
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